General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



DQE/NASA/0207-79/1 
NASA TM-79052 

(NASA-TN-790S2) PLOW FRICTION OP THP 
T'JRP'ITPNT COOLANT PLOW IN CRYOGENIC POROTS 
CABLPS (NAPA) Ifl p HC A02/«F A01 C3CL 20D 

G 3/34 

FLOW FRICTION OF THE 
TURBULENT COOLANT FLOW 
IN CRYOGENIC POROUS CABLES 


N79-20341 

One Las 
17292 


Robert C. Hendricks 

National Aeronautics and Space Administration 
Lewis Research Center 

and 

V. M. Yeroshenko, L. I. Zaichik, and L. S. Yanovsky 
Krzhizhanovsky Power Engineering Institute 


Work performed for 

U. S. DEPARTMENT OF ENERGY 
Office of Energy Technology 
Division of Electric Energy Systems 



TECHNICAL PAPER 

Fifteenth' International Congress on Refrigeration 
sponsored by the International Institute of Refrigeration 
Venice, Italy, September 23-29, 1979 



DOE /NAS A/0207 -79/1 
NASA TM-79052 


FLOW FRICTION OF THE 
TURBULENT COOLANT CABLES 
IN CRYOGENIC POROUS CABLES 


Robert C. Hendricks 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

and 

V. M, Yeroshenko, L L Zaichik, and L S. Yanovsky 
Krzhizhanovsky Power Engineering Institute 
Moscow, U. S. S.R. 


Prepared for 

IIS. DEPARTMENT OF ENERGY 
Office of Energy Technology 
Division of Electric Energy Systems 
Washington, D.C. 20545 

Under U. S. - U. S. 5. R. Energy Agreement Subtask 080207 


Fifteenth International Congress on Refrigeration 
sponsored by the Inter national Institute of Refrigeration 
Venice, Italy, September 23-29, 1979 


6986 


PAGE 1 


FLOW FRICTION OF THE TURBULENT COOLANT FLOW 
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Cleveland, Ohio 44135 

and 

V. M. Yeroshenko, L- I. Zaichik, and L. S. Yanovsky 
The Krzhizhanovsky Power Engineering Institute, 

Moscow, U.S.S.R. 


ABSTRACT 

Considered are cryogenic power transmission cables with 
porous cores. Calculations of the turbulent coolant flow 
with injection or suction through the porous wall are 
presented within the framework of a two-layer model. 
Universal velocity profiles have been obtained for the 
viscous sublayer and flow core. Integrating the velocity 
profile, the law of flow friction in the pipe with injection 
has been derived for the case when there is a tangential 
injection velocity component. The effect of tangential 
velocity on the relative law of flow friction is analysed. 
The applicability of the Prandtl model to the problem under 
study is discussed. It is shown that the error due to the 
acceptance of the model increases with the injection 
parameter and at lower Reynolds numbers; under these 
circumstances, the influence of convective terms in the 
turbulent energy equation on the mechanism of turbulent 
transport should be taken into account. 
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INTRODUCTION 

One of the best techniques foe cooling cryogenic pouer 
transmission lines consists in the utilization of cores 
permeable for the coolant. When cables of this type are 
involved, the turbulent coolant flow in a long pipe with 
injection or suction of the coolant through the porous wall 
should be considered in order to predict thermal ?_"d 
hydraulic behaviors. 

Previous experimental and theoretical investigations 
have been mainly focused on boundary layer flows while flows 
with injection or suction have received inadequate attention 
in literature. It is a common practice to estimate 
parameters of turbulent flows in channels with injection or 
suction using procedures developed for the boundary layer 
(for instance, ref. 1). Besides, numerical calculations for 
flows with injection or suction have been reported (refs. 2 
and 3) . 

However, for practical purposes, it is of great 
interest to develop calculation procedures applicable to 
channels with injection or suction and involving relatively 
simple analytical equations. Such a procedure is proposed 
in this paper. 

The flow of an incompressible fluid with constant 
physical properties through a circular pipe is described by 
the following equation 


3 2,3 

ru 


1 9P _ 1 3_ 
p 3x p 3r 


( 1 ) 


When normal and tangential injection velocities are 
taken into account, the integration of equation (1) over the 
pipe cross-section yields 


4- BU 2 + — ^ i dP _ £l2Z. 

dx R p dx pR 


( 2 ) 


where U s (2/R 2 )/ urdr 
0 


velocity, f3 £. (2/R 2 U V u 2 rdr 

0 


is the mean mass 


is the coefficient of impulse flow. - 

Since within the turbulent flow core the velocity 
profile is nearly uniform, it may be approximately assumed 
that for the entire cross-section of the pipe exclusive of 
the region near the wall 
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Using equations (2) ana (3) , equation (1) can be rewritten 
as follows 
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(4) 


Below, in accordance with the assumption of local 
similarity (ref. 2) , u+ and t + are presumed to be 
dependent only upon the transverse coordinate r. 

Moreover, bearing in mind experimental results obtained 
with injection and given in reference 4 ds well as numerical 
calculations carried out for suction in reference 2, we 
assume the linear dependence between the radial velocity and: 


v “ v r 
w 


(4a) 


Now the solution of equation (4) is 



To describe the turbulent impulse transport, the Prandtl 
model may be employed 



where l - Kyf(y), f(y) 1 at y + 0 

Considering equation (5) , the integration of equation 
(6) gives the velocity profile 



-± 

k J yf (y) 


dy 


(7) 


As in reference 5, y Q is a distance from the wall to where 
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the velocity equals to u . _ 

Subtracting from equation (7) its expression at y = 1, 
one obtains 



F(y) 


( 8 ) 


The above equation gives the velocity defect for pipes 
with in jection/suction. _ 

According to reference 5, the explicit _ function _F(y) 
may be given as follows 


F(y) - A = i - 2.5 


(9) 


From equations (8) and (9) , it may be concluded that 
for the entire pipe cross-section the following logarithmic 
dependence is valid and equation (7) becomes 



n 



Alny + + B 


( 10 ) 


where B = 5.5 - constant independent of injection. 

With normal injection, equation {10) for the region 
near the wall is well known (cf. refs. 1 and 6) while for 
the flow core it can ba easily found from equations {8) and 
(9). 

The explicit expression for the velocity profile can be 
derived from equation (10) 


u, 


u + Alny + + B + ~ (Alny + + B) 2 
+ 


( 11 ) 
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The presence of a tangential velocity, the velocity 
profile vithiu the viscous sublayer of turbulent flow with 
injection is 


u . " u + 
+ 


ex P(V+) - 1 


V, 


( 12 ) 


Within the framework of two-layer flow model, the most 
natural manner for determining the viscous sublayer 
thickness consists in setting the right-hand sides of both 
velocity expressions (eqs. {11) and (12)) equal to each 
other 


+ 2 
Aln6 + + B + ~ (AlnS + + B) 


exp(V + 6 + ) - 1 


V. 


(13) 


Equation (13) shows that the thickness 6 + is 
independent of the tangential velocity u w+ . 

Figure 1 presents the calculated dependence between 6 + 
and V, . Moreover, the experimental data from reference 8 
along + with the viscous sublayer thickness obtained in 
references 9 and 10 are also shown in figure 1. 


«+ ‘ v7 M 1 + °o v +) 

*T* 


(14) 


where a ^ = 13.1 according to reference 9 and c^ =11.9 
a ccor d ing^ to reference 10; V + is greater than '-a” 1 • — 

It should be noted that a measure of stability, Reg 
including parameters calculated from equations (12) and (13) 
for the viscous sublayer boundary, reduces as the injection 
velocity grows and it increases with the suction velocity. 
This is in agreement with the physical picture of the 
injection influence on the region near the wall (refs. 8 and 
10 ). 

Employing the relation u = /8/A , the friction 
resistance in the pipe can be + derived. U + is given by 
equation (12) for the viscous sublayer and by equation (11) 
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for the remainder of the flow. 

The friction resistance in a pipe with slanted 
injection (suction) can be found by integrating the veLocity 
profile to give: 




e*p(V + ) ^P(V+) 

«2_2 + V.R, 




(15) 


where 



The terms of an Qrder_.O.f . <S^/R + are omitted in 
equation (15) since 5 + /r + « 1 „ ' 

In the absence of x n jection or suction, equation (15) 
gives the flow friction for a round tube 




“ AlnR 
X 0 *0 


+B ’! A+ rTC 


6 , + 2A - 2Aln6 

+ 0 + 0 


- 2B 
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(16) 
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where 


R 


Re 


■ + o 2 


*/s 

f 8 * 


6 = 11,64 

*0 


Solving equation (16) by the method of intention at 
high Re numbers, the explicit expression for A(Re) can be 
obtained. 

Four interations produce 
- (l + ||) /&ln I s - + D - Ain j^Aln |- + D - Ain [Ain + D - Ain (Ain | 


where 


D - B ~ 


(17) 


N - 6 (6 + 

*0\ 0 


OV ’0 
3A 


2A - 2Aln5 - 2B 

< A 


It is convenient to present calculations as a dependence 


between the relative flow friction 


X/X 


.0 


and injection or 


suction parameter b « 8v /UXq which is a common practice in 
the boundary layer theory (ref. 1). 

When Re tends to infinity, it follows from equations 
(15) and (16) that 


11m 
Re ■+ “ 


k ■■) ■ <■ ■ ••)’(* - tj 


(18) 


where 


% * -H 1 - “«)• “w - TT 


absence of the 
(10) coincides with 
by Kutateladze- 


It can be readily seen that in the 
tangential injection velocity, equation 
the limit friction dependence suggested 
Leontyev for turbulent boundary layers. 

Figure 2 shows ’J'(b) calculated from equations (15) 
and (17). 

The contribution of the viscous sublayer to the flow 

friction appears to be insignificant even at Re = 10^ and 

decreases further as the Reynolds number grows. 

When there is an injected flow in the pipe, the flow 

friction decreases with Re, similarly to the case with the 

boundary layer (ref, 1) , while in the presence of suction 

flow, the Reynolds number influence becomes insignificant. 

The latter is in agreement with reference 2 and the 

simplified equation of Wallis (ref. 2-Discussion) . 

Figure 2 indicates that the results obtained are in 

satisfactory agreement with other experiments and 

calculations carried out for injection and suction 

conditions at m - 0 . 

w 


it 


Figure 3 'illustrates the influence of the tangential 
velocity on the relative flow friction in. a pipe with 
injection or suction.. The skin friction stress Y 
decreases smoothly as the tangential velocity component u w , 
directed along the flow, increases; while Y increases when u 
directed against the flow, increases. 

The procedure proposed hero is based on certain 
assumptions the main of which is the acceptance of the 
Prandtl model (equation (6)). To estimate the region 
wherein equation (6) holds, it is necessary to consider the 
equation of turbulence energy balance without the diffusion 
term but containing terms characterizing couvection, 
generation, and dissipation (ref. 15) 


u 




l „ au 
p 3r 



(19) 


where i ~ turbulence scale assumed to be equal to the 
Nikuradse mixing region length 


_ — 2—4 

i = 0.14 - 0.08r - 0.06r (20) 

Similarly to the assumption of local similarity made in 
above in respect to t, and u + , the turbulent energy, e , 
is assumed to be a function only of the transverse 
coordinate, E + » E + (r) * 

To determine the Reynolds stresses, the Kolmogorov 
correlation can be used (ref. 15) 


T 


+ 


a£E 


1/2 



( 21 ) 


If the convective terms are absent in equation (19), 
the turbulent tangential stress and energy can be obtained 
from equations (19) to (21) 



( 23 ) 
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The constants o and c may be determined from 
equations (22) and (6), a/a/c » 1 t and, since in the region 
near the wall, the ratio between the turbulent tangential 
stress and turbulent energy remains constant, this 
pt-rarae tar, according to reference 16, can be taken to be 
t,/E » (ac)l/ 2 =* 0.26 T 
^ Consequently, a » 0.51, c « 0.13 . 

The assumption of turbulent energy local similarity makes it 
possible to define the term u 8E/3x in equation (19) 


9_E 

3x 



uE 


+ dx 


2uu * E + 



or, taking into account the continuity equation, one obtains 


u 


3E 

3x 


n 


R 




(24) 


Using eguations (4a), (21), and (24), equation (19) can 
be written as follows 



-2V, 




pi/ 2 
'+ T 


- dE + 

^ dT 


, 1 / 2 “ 


(25) 


where the velocity profile u + is determined from equation 
(11) . The iteration procedure lias been employed to solve 
equation (25). Considering the right-hand side as a 
disturbing factor, the solution of equation (25) in the 
first and second approximations may be written 



Tho approximation E+^ is justified under the 
assumption that tho turbulence onorgy generation and 
dissipation are equal, that is, when tho Prandtl model is 
valid. , 

Calculated values of E 1 and E 2 at Re «* 10 and 
w w = 0 normed on u^ 0 are plotted In figure 4. 

It is evident that the maximum of turbulence increases 
and moves away from the wall as the injection flow grows. 
This is found to be in qualitative agreement with 
experiments reported in reference 17 for a rectangular duct 
with the injection flow entering the duct from one of the 
vails. At high parameters of suction, the flow in the pipe 
center tends to become laminar. 

Discrepancy between E 1 and E., due to the influence 
of convective terms on the turbulent transport mechanism, 
characterizes the region of applicability of results on skin 
friction in pipes with injection or suction. 

As the injection parameter increases (b < -4) and Re 
decreases, the error due to the Prandtl model grows and one 
should consider the influence of convective terms in the 
turbulent energy equation on the turbulent transport 
mechanism. 


NOMENCLATURE 


x,r 

p 

p 

t 

R 

5 

v 

>j* - A w /p 
A 

Dimensionless 


axial and radial coordinates 

axial and radial velocity corresponding 

pressure 

density 

shearing stress 
pipe radius 

viscous sublayer thickness 
kinematic viscosity 
shear velocity 
skin-friction coefficient 
variables: 


by 


y “ 1 - r 
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y u * 



Re 


?RU 

v 


R. 




Subscripts: 

W at the wall 

0 without injection or suction 


ACKNOWLEDGEMENT 


The assistance of J. A. Hendricks and J. R 
preparing this document is greatly appreciated. 


Hendricks in 





PAGE 12 


REFERENCES 

1. Kutateladze, S. S. , and A. I. Leont'yev, Toplomassoobmen i 

treniye v turbulentnom pogranichnom sloye , (Heat and 
Mass Transfer and Friction in a Turbulent Boundary Layer) , 
Moscow, Energiya Prose, 1972. Also NA3A-TTF-805 , 1974. 

2. Kinney, R. B., and Sparrow, E. M. , ’Turbulent Flow, Heat 

Transfer and Mass Transfer in a Tube with Surface Suction, 1 
ASME Trans., Heat Transfer , Ser. C, Vol, 92 (1970), pp . 
117-125 . 

3. Sviridenkov, A. A. and lagodkin, V. T. , 'Flows in Initial Re- 

gions of Channels with Porous Walls,’ Izvestiya Akademli 
Nauk SSR, Mekhanika zhidkosti i gaza , Mo. 5 (1976), pp. 

43-48. Translated Fluid Dynamics , Vol. 11 (1976), pp. 
689-693. 

4. Olson, R. M., and Eckert, E. R. G., 'Experimental Studies 

of Turbulent Flow in a Porous circular Tube -ith Uniform 
Fluid Injection through the Tube Wall,' ASME Trans., 

Applied Mechanics , Ser, E, Vol. 33 (1966), pp. 7-17. 

5. Schlichting, H. , (Boundary Layer Theory) , Sixth ed. , McGraw- 

Hill, New York (1969). 

6. Dorrance, W. H., and Dore , F. J. , 'The Effect of Mass Trans- 

fer on the Compressible Turbulent Boundary Layer Skin 
Friction and Heat Transfer,' J. Aeronautical Sciences, 
ol. 21 (1954), pp. 404-410. 

7. Kutateladze, S, S., Leont’ev, A. I., and Mironov, B. P., 

Turbulent Boundary Layer on a Porous Wall with a Differ- 
ential Pressure Gradient in the Region of Strong Air 
Blasts and Finite Reynolds Numbers,’ Teplo- Massoperenos , 
Dokl. Vses. Soveshch., 3rd., (Heat and Mass Transfer) , 

Vol. 1 , Lykov, A. V., ed. , Moscow, Energiya Press, (1968), 
pp . 24-32. 

8. Ermakov, A, L., Eroshenko, V, M. , Terent'yev, Yu. N., and 

Yanovskli, L. S., 'Investigation of Laminar Sublayer on 
Permeable Surfaces with Injection,’ Izvestiya Akademii Nauk 
SSSR, Mekhanika Zhidkosti i Gaza , No. 5 (1977), pp . 66-72. 

English Translation, Fluid Dynamics , Vol. 12 (1977), 
pp. 696-702. 

9. Rubesin, M. W. , 'An Analytical Estimation of the Effect of 

Transpiration Cooling of the Heat-Transfer and Skin-Friction 
Characteristics of a Compressible Turbulent Boundary Layer, 1 
NACA Technical Note 3341 (1954). 

10. Motulevish, V. P., Turbulent Heat and Mass Transfer on a Plate 
During Porous Drawing Down and Injection of Different Gases,’ 
Inzhenerno-Fizicheskiy Zhurnal , Vol. 6, (1963), No. 1, pp. 

3-13. English translation, J. Engineering Physics , Vol. 6 
(1963) . 


.W~U 




11. Kochenov, X. S. and Novos,. 3 * skii , 0. Yu., Hydraulic Re- 

sistance of Channels with Permeable Walls, Inzhenerno- 
Fieicheskiy Zhurnal , Vol. 16, Ho. 3, 1969. pp. 405-412; 
Translated J. Engineering Physics , Vol. 16 (1969), pp. 
275-281. 

12. Elana, M., 'Effect of Drawing Down on the Average Velocity 

and Temperature During Turbulent Flow in a Round Tube,' 

In the collection; Teplo- i massoperenos (Heat and Mass 
Transfer), Vol. 9, Part 2, Minsk, 1972. 

13. Burnage, H., 'Contribution A L'Etude De L ' Ecoulemen t Turbu- 

lent Dans Une Conduite Cylindrique A Paroi Poreuse' (Con- 
tribution to the Study of Turbulent Flow in a Cylindrical 
Pipe With Porous Walls), Paris, Service de Documentation 
et d ' Inf ormation Technique de 1 1 Aeronautique , Publication 
Scientifique et Technique N. T. 114, (Aug. 1962). 

14. Wallis, G. B., 'Pressure Drop with Lateral Mass Extraction. 

I. Pressure Gradients for Air Flowing Along Porous Tubes 
with Uniform Extraction of the Walls,' Proc. Inst. Mech- 
anical, Engineers (London) , Vol. 180, (1965-66), pp. 27-35. 

15. Monin, A. S, and Yaglom, A. M. , Statisticheskaya gidromek- 

hanika (Statistical Hydrodynamics) , Vol. 1, Moscow, Nauka 
Press, (1365). Also JPRS-37763. 

16. Laufer, D. , 'The Structure of Turbulence in Fully Developed 

Pipe Flow,' NACA Technical Report 1174, 1954. 

17. Paleev, I. I. Agafonov, F. A., and Dymant, L. N., 'Experi- 

mental Study of Isothermal Turbulent Flow in a Rectangular 
Channel with a One-sided Air Injection,' Izvestiya vysshikh 
uchebnykh zavedeniy SSR, seriya Energetika, No. 1, (1970), 

pp . 65-70. 


FRICTION DE L ' ECOULEMENT TURBULENT EN CABLE POREUX CRYOGENIQUE 

SUMMAIRE; Ici sont considerais les cables de puissance 
cryogenique avec injection ou aspiration a travers la muraille 
poreuse, est presents dans le cas d'un modele avec deux couches. 
Le porfile de velocite universelle est derive pour la couche 
vi?'.ause et le centre d ' ecoulement . Avec 1' injection du 
profile de velocite, ont arrive a la loi qui directe la friction 
de 1' ecoulement dans un tube avec injection pour le cas ou il 
y a un component de velocite d' injection tangentiel. L'effet de 
la velocite. tangential sur la loi relative a 1'coulement de 
friction est analyse. L 1 application du modele de Prandtl au 
probleme ici est discoute. II est montre gue l'erreur du modele 
est augmente par le paramete d' injection pour les numbres 
Reynolds plus bas ; sons ces circumstances, 1' influence des 
termes convectives dans 1' equation d'energie turbulente devra.it 
etre pris en comptes. 
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Figure 1. - Viscous sublayer tnicltness vetsus in- 
jection or suction velocity. 
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tive flow friction. 
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Figure 3. - Injection or suction tangential velocity versus flow friction. 
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Figure 4. - Distribution of turbulent energy in the pipe 
cross-section. 



